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Abstract 

A joint characterisation of the observability and controllability of a 
particular kind of discrete system has been developed. The key idea of 
the procedure can be reduced to a correct choice of the sampling sequence. 
This freedom, owing to the arbitrary choice of the sampling instants, is 
used to improve the sensitivity of system observability and controllability, 
by exploiting an adequate geometric structure. Some qualitative examples 
are presented for illustrative purposes. 

1 Introduction 

The concepts of observability and controllability, introduced first by Kalman 
[B], play an important role in modern control theory. In fact, these properties 
often govern the existence of a solution to an optimal control problem. 

Observability and Controllability of discrete-time systems have also been 
treated in the literature in a generalised form. A survey of the main results for 
discrete system sampled nonuniformly can be found in |7j. 

In general, both concepts have been characterised by criteria mutually inde- 
pendent, but, in this work, it is demonstrated that, if certain conditions not very 
restrictive are imposed on the continuous system, then the characterisations of 
observability and controllability for the discrete system can be unified. 

There is an important problem which arises in discrete-time systems but 
not in continuous-time systems. A linear system which is completely observable 
and controllable may lose these properties after the introduction of sampling. 
This stresses the importance of the sampling sequence to guarantee the above- 
mentioned internal properties. 
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There are restrictions on the aperiodic samphng sequence, stronger than in 
the periodic case. At any rate, this partial freedom in the choice of the samphng 
instants can be conveniently used to obtain a rather well conditioned system, a 
reduction in the propagation of measuring and/or rounding errors etc. 

Further interesting applications of these ideas, e.g. the sensitivity of system 
observability and controllability, are treated by exploiting an adequate geometric 
structure. In this way, the familiar analytical techniques [7] can be presented in 
a more intuitive form. 

2 Basic assumptions and general considerations 

This discussion is restricted to: 

a) linear time-invariant single- input / single-output differential systems of finite 

order n 

b) continuous-time systems completely controllable and observable in the sense 

given by Kalman [S] 

c) systems whose transfer function is a strictly proper rational function. 

Consequently, their impulse response h(t) will be a particular solution of an 
nth-order homogeneous linear differential equation with constant coefficients of 
the form 

h^"\t) + aih^^'-^^t) + . . . + ar,h{t) = t>0 (1) 

Therefore 

n 

= ^^(^,(0 (2) 

1=1 

where Ci G C are constant coefficients. ipi;R C (i = 1, . . . , n) is the 
fundamental system of solutions of eqn. [l] 

In the state space, the systems considered can be described by the following 
equations 

X{t) = AX{t) + bu{t) [Xq = X{Q)) 

Y{t) = cX{t) (3) 

where X e M" denotes the state- vector and u,y & M. are the scalar input 
and output, respectively. The matrices A,b,c are of appropriate orders and 
constant. 

Canonical realisations in the state space obtained from the impulse response 
will be used. In particular, 

(i) Observability canonical form {Aob,bob, Cob)- 
Aob is a n X n bottom-companion matrix with 

(-a„, -a„_i, . . . , -ai) (4) 



2 



in the last row, 

bob = {hi,h2, . ■ . , hn)' 
(' denotes the transpose) where 



hi+i — 



d^h{t) 



(z = 0,...,n-l) 



(5) 



(6) 



correspond to the n-first Markov parameters of the impulse response h(t) 



and 



Cob = (1,0,...,0) 
(ii) Controllability canonical form {Aco^b^o-^Cc 



(7) 



A' 



bro = C' 



ob 



(8) 
(9) 

^ob (10) 

All minimal realisations are related by similarity transformations and, in 
each particular problem, the most adequate one will be selected. 



= b', 



3 Joint characterisation of observability and con- 
trollability of nonuniformly sampled discrete 
systems 

The observability (controllability) of a discrete system depends on the observ- 
ability (controllability) of the continuous-time system, plus some additional con- 
ditions on the sampling sequence. 

The problem of observing (controlling) the state of any realisation by means 
of the sampling can be reduced to the task of solving a system of linear equa- 
tions. Consequently, an adequate choice of the sampling instants guarantees the 
compatibility of this system. 

The question of the controllability will be discussed first. The corresponding 
results for the observability will be given later. Finally, a joint characterisation 
of both internal properties will be presented. 

3.1 The controllability problem 

Let ( A, b, c) be an arbitrary minimal realisation of order n for the kind of sys- 
tems under study. The system considered will be completely n-controUable 
(controllable in n steps) [T] , [7] if, for any initial state Xq of eqn. |3] the system 
can be directed to X = by means of n impulse inputs applied at n consecutive 
sampling instants. The solution of the state-space equation for the system eqn. 
[3] at time t„ can be written as 

n-l 

X{tn) = eXp{Atn)Xo + ^ GiUi 
i=0 



3 



exp{Atr,)Xo + [Gn- 



J Go] 



uo 



(11) 



where 



G, = exp{A(t„-t,+i)) exp{A(T,+i~r))xbS{r~t,)dr {i = 0, . . . ,n-l) 

J-o 

(12) 

Ti+i — ti^i — ti (13) 

is the length of the samphng interval between two consecutive sampling 
instants. The scalar input is 



u{t) — 6{t — ti)ui (i = 0, . . . , n — 1) 



(14) 



Ui being the value of the impulse input at the sampling instant ti . To consider 
the transference from an initial state to a final state in n steps, we study the 
rank of the matrix [G„_i, . . . , Go]. Indeed, 

Gi = B exp{J {tn - tn^i)) exp( J (t„_i - ti))yo (15) 

where J is the Jordan canonical form of matrix A and B is the invertible 
matrix of the change of basis: 



Therefore, we must compute the value of 

Det = [yo,exp{J{tn-i - i„-2))2/o, ■ ■ ■,exp{J{tn-i - ta))yo] 
or briefly. 



with 



Det ^ [exp{J{ajn)yo] (m = 0, . . . , n - 1) 

(ao = 0) 



Q^m ^n — 1 ^n—ni—l 

By similarity transformations 

yo = B-'b = B-^X 



(16) 

(17) 

(18) 

(19) 
(20) 



where B^o is the matrix of the change of basis of A^o to the Jordan canonical 
form. Then we denote the components of ya by 



/ I 1 r r \/ 

yo — [Vl 1 ■ ■ ■ 1 J/ml J • ■ • I 2/l J • • ■ I ymr) 



(21) 



mj [j = 1, 



") being the multiplicity of the r different eigenvalues of 



matrix A with r < n. It is easy to see [5] that, for the controllability canonical 
form, 
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Bpgbco — (Ci, C2 



,Cn)' 



(22) 



with [Ci) defined as in eqn. 2; thus making use of the Laplace's expansion 
by minors, the determinant eqn. 17 can be factorised as follows: 



Det[exp{Jam)yo] = N1N2 Det[ipi{am)] 



(23) 



1 
0! 



1 

(mi - 1)! 



1 
0! 



1 

{mr - 1)! 



(24) 



Ni = Del 



Vml 



Note that N2 will be non null if, and only if, 



(25) 



J/L^O (i = l,...,r). 



(26) 



which is guaranteed, according to the previous significance of the components 
of t/O) because only minimal realisations are considered. 

Finally, [(/3i(am)] (i = 1, • • . , n; m = 0, . . . , n — 1)] is an n x n matrix of the 
form 



</5i(ao) 
¥'1(0:1) 

(Pi(q!„_i) 



(27) 



{ipi{a)) being the fundamental system of solutions of eqn. 1. The value 
of Det[(pi{am)] will depend on the choice of the sampling instants. Let us 
now consider one special aspect here, because it can be the source of some 
terminological problems. 

The nonsingularity of the matrix [G„_i, . . . , Go] is a necessary and sufficient 
condition to ensure that any initial state Xq can be taken to an arbitrary state 
Xn in n steps: this is called n-reachability. If the state X„ coincides with the 
origin, then the above condition is sufficient to insure that any initial state Xq 
can be taken to the origin in n steps: this is called n-controllability. 

The n-reachability is a more restrictive condition than the n-controllability. 
In fact, the n-reachability implies n-controllability, but the converse may not 
be true. 
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From eqn. 23 and the previous considerations, the following result is derived: 
Lemma 1: The realisation (A, &, c) is completely n-reachable (reachable in n 

steps) if, and only if, n consecutive sampling instants are chosen in such a way 

that 

Det[ip,{ara)] ^ Q (i = l,...,n; m = 0,...,n-l) (28) 

Proof: The proof is evident from the factorisation of the expression 23. 

Note that, for this kind of system, the n-reachability depends on the char- 
acteristic modes of the continuous system and on the choice of the sampling 
instants. 

We remark that the hypothesis (h) in Section 2 is necessary, because, if 
this condition is not verified, then the canonical realizations obtained from the 
impulse response would not be actual realisations in the state space for the 
system considered. Consequently, the above result would not be true. 

In eqn. 14, the scalar input is defined as impulse inputs at the sampling 
instants. In practice, a control of the form 

u{t) = u{ti) ^Ui U<t < U+i (29) 

is generally used. 

For technical reasons, u{t) must be generated with the aid of a filter, so that 
the equations of the device which generates u{t) can be included in eqn. 1. The 
results are completely analogous. 

3.2 The observability problem 

The observability problem can be discussed without loss of generality putting 
u(t) = 0. For the same arbitrary minimal realisation (A, h, c) the observability 
question is studied. The system considered will be completely n-observable 
(observable in n steps) [T], [7], if any initial state of eqn. 3 can be calculated 
from n values of the output taken at n consecutive sampling instants. 

Therefore, according to eqn. 3 and for the same values ao, ai, . . . , ctn-i as 
before, the following system of linear equations is set: 

v{am) = cexp{Aa„i)Xa (m = 0, . . . , n - 1) (30) 
In mathematical terms, the condition of n-observability means that 

rank[cexp{Aam)] — n (m = 0, . . . , n — 1) (31) 
The linear system eqn. 30 can be rewritten as 

yiam) ^ cB exp{Jam)zo (m = 0, . . . , n - 1) (32) 

where J and B are defined as before 

zo = B-'Xo (33) 



6 



Now, we study the value of 



Det[cB exp{Ja„i)] (m = 0, 
By similarity transformations 

cB = Cob Bob 



1) 



(34) 



(35) 



where Bob is the matrix of the change of basis of Aob to the Jordan canonical 
form. It is known j5j that, for the observability canonical form, 



(36) 



mi 

CobBob - (1,0^-, 0, 1,0,..., d) 

, r) is the multiplicity of the eigenvalues of matrix A, 



where nij {j = 1, 
with r < n. 

Thus, making use of the Laplace's expansion by minors, the determinant 
eqn. 14 can be factorised as follows: 



Det[cB exp{Jam)] = M1M2 Det[ip,{am)] 



where 



^ 0! 



1 



(mi-1)! 
h 



1 
0! 



{rrir - 1)! 



Mo = Det 



h 



Ir 



(37) 
(38) 

(39) 



with Ij (j : 



l,...,r) identity matrix of dimension rrij x nij. Finally, 
■ ,n), (m = 0, . . . , n — 1) is the same matrix as the one defined 
in eqn. 27. Now, from the expression 37, the following result is derived: 

Lemma 2: The realisation {A, b, c) is completely n-observable (observable in 
n steps) if, and only if, n consecutive sampling instants are chosen in such a 
way that 



Det[ipi{arn)] 7^ {i = l,...,n; m = 0, ...,n-l) 



(40) 



Proof: The proof is evident from the factorisation of the expression 37. 

Note that, for this kind of system, the n-observability depends on the char- 
acteristic modes of the continuous system and on the choice of the sampling 
instants. We remark that the result is analogous to that of the preceding sub- 
section. 
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3.3 Main result 



The results obtained in the two preceding subsections can be summarised as 
follows: 

Theorem: A system verifying the conditions (a), (h) and (c) \n Section 2, 
is jointly n-reachable and n-observable, if, and only if, n consecutive sampling 
instants are chosen in such a way that 

Det[ip^{a„i)\ ^ Q (i = 1, . . . ,n; m = 0, . . . ,n- 1) (41) 

Proof: The proof is evident from the use of the results obtained in Sections 
3.1 and 3.2. 

The condition 41 imposes a rather weak restriction for the choice of the 
sampling instants. In fact, intervals can be specified [3], [7] so that complete 
reachability (observability) is preserved. 

We can also remark that, for this kind of system, the n-reachability and 
n-observability are inseparable concepts. These systems are either n-reachable 
and n-observable or, if not, they are neither n-reachable nor n-observable. 

It must be noticed that the condition 41 depends exclusively on the charac- 
teristic modes and the sampling sequence. Thus, the n-reachability (n-observability) 
for these systems is independent of the chosen realisation. 

If the system does not verify the condition (h ), then the results obtained will 
still be valid for the subsystem controllable and observable. 

Note how reachability involves statements about inputs and states, while 
observability about outputs and states. Nevertheless, both aspects are related 
to the same intrinsic system properties. Consequently, the algebraic character- 
ization for reachability and observability can be unified. 

4 Problems of sensitivity of system observabil- 
ity and controllability 

In stability theory, it is not enough to know whether a system is stable or 
unstable but also the degree of stability. In a similar way, it is important to 
know the degree of observability (controllability) of the system; in the sense 
that there are different levels of certitude in the process of resolution of the 
corresponding linear equations. 

From a heuristic viewpoint, the degree of observability (controllability) de- 
pends on the spatial relationship among the column vectors of the observability 
(controllability) matrix. In particular, maximum observability (controllability) 
is obtained when these vectors are mutually orthogonal [7]. 

Different geometric structures, which show the evolution in time of the above 
mentioned vectors, can be found in [3]. Thus, the question is to choose conve- 
niently the sampling instants to obtain maximum orthogonality. 

We start from an arbitrary minimal realisation {A, b, c) of the given system 
eqn. 3. The controllability matrix 
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Figure 1: Relationships among the vectors Yo,Yi 



[Gn-i,...,Go] (42) 
can be easily reduced to the matrix 

[Fo,Fi,...,F„_i] (43) 



Yi = exp{J{tn-i-tn-i-i))yo (z = 0,...,n-l) (44) 

( J, yo defined as before) more convenient for the geometric interpretation; as 
the factor B exp{J{tn — tn-i)) in eqn. 15 affects only the module of the column 
vectors but not their spatial position. Analogous results can be obtained making 
use of the transpose observability matrix. 

Note that matrix 43 depends exclusively on the system characteristic modes 
and the sampling sequence; thus the degree of observability and controllability 
is independent of the chosen realisation. Now the problem is reduced to a right 
choice of the sampling instants, in such a way that the vectors Yi are mutually 
orthogonal. Some qualitative examples are presented: 

Example 1: We are going to consider a 2nd-order model with a pair of 
complex eigenvalues 

a + i6 e C (6 > 0) (45) 

then 

\Yq,Yi] = [exp{Jao)yo,exp{Jai)yo] (46) 



a —b 
b a 



(6>0) 



(47) 



ai,yo defined as before. 
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As we are in R^, the geometric interpretation is very simple. The generic 
operator exp{Ja) apphed to the vector yo can be viewed as follows. It is a coun- 
terclockwise rotation through ba radians, followed by a stretching (or shrinking) 
of the length of t/o by a factor exp{aa) [1] . 

From this interpretation, it is easy to see (Fig. 1) that the vectors Yq,Yi, 
will be mutually orthogonal if, and only if, 

5ai = 6(ii-to) = (2m + l)| (m = 0,l,...) (48) 

In this case, maximum observability (controllability) is achieved. In this 
situation, expression 48 coincides with the results obtained by Troch [7| for the 
same model, when the problem of a minimum transmission of the measuring 
errors is discussed. 

Example 2: Let us now study a 3-order model with a real pole and a complex 
pair: 

A e M, a + jbeC {b>0) (49) 

Therefore, then 

[Yo,Yi,Y2] = [exp{Jao)yo, exp{Jai)yQ, exp{Ja2)yQ] (50) 

with 



J 



a -b 
b a 
A 



[b > 0) (51) 



ai , yo as usual. The problem is treated in the 3-dimensional space and the 
geometric structure can be described as follows: the generic vector Y{a) is 
written as 

Y(a) = {exp{aa) cos{ba) , exp{aa) sm{ba) , exp{Xa)y (52) 
thus it is a spiral on the surface revolution: 

Z^{X^+ y2)V2a (53) 

whose form is determined by the real part of the system eigenvalues. 

The vectors Yq, Yi,Y2 have their origin at the point (0, 0, 0) and their end at 
the points Y{ao),Y{ai),Y{a2) of the parametric curve (Fig. 2). Heuristically 
we can imagine the same vectors as before (2-dimensional example) pointing 
upwards from the XY-plane as they have a third component on the Z-axis. 

Given Yq and Yi , the question is now to choose the sampling instant such that 
the vector Y2 obeys its motion law and, simultaneously, is the most orthogonal 
vector to the other two. 

In fact, let Po,Pi,F2 be the projection on the XF-plane of the vectors 
Fq , ^1 , ^0 X ^1 7 respectively: 

Po = {xq, yo) = {exp{aao) cos{bao) , exp{aao) sm{bao)) (54) 
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Figure 2: Relationships among the vectors Fg, Yi and Y2 



Pi = {xi,yi) = {exp{aai) cos{bai) , exp{aai) sui{bai)) (55) 
P2 = {X2,y2) = {yozi - yiZo,xiZo - xqZi) (56) 

with 

Zi = expiXai) {i = 0, 1) (57) 
The sampHng sequence must be selected such that 

ba2 = arccos ^ ° ' ^ ^ = M (58) 
I "o| I "2I 

then 

b{t2 -to)=M + 27rm (m = 0, 1, . . .) (59) 

This number can be understood as the angle of a comiterclockwise rotation. 
To obtain the more adequate value for m, we analyse the geometry of this 
situation. Let Q2 be a vector orthogonal to Yq and Yl with its end on the 
revolution surface. Indeed, 

Q2 = fi{Yo X Yl) (60) 
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with /X easily computable according to eqn. 53; then the best value of m will 
be the integer for which the expression 



{li\xl + _ exp{^{M + 27rm)) 



(61) 



is minimum. The process is computed again for each new sampling instant. 

These examples show that the optimal sampling sequences involve only dif- 
ferences between sampling instants, but not absolute positions of such instants 
on the time axis. It perfectly agrees with the kind of time invariant systems we 
are considering. 



5 Conclusions 

It has been shown that, for completely observable and controllable systems 
described by eqns. 3, a joint characterisation of the n-reachability and n- 
observability of the corresponding discrete systems can be given. The formula- 
tion considered stresses the importance of the sampling instants to preserve the 
aforementioned internal properties. In this way, such systems have an additional 
element for analysis and manipulation. 

There are no overrestrictivc conditions on the choice of the sampling se- 
quence. Moreover, it has also been possible to give strategies for some special 
cases. At the same time, different geometric structures have been used to im- 
prove the sensitivity of system observability and controllability, according to a 
correct selection of the sampling instants. These geometric structures can also 
be used with other performance criterions to give some insight on the system 
evolution in terms of familiar concepts. 
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